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We compare the BCS and Overhauser effect as competing mechanisms for the destabilization of 
the quark Fermi surface at asymptotically large chemical potential, for the special case of 2 space 
and 1 time dimension. We use the framework of perturbative one-gluon exchange, which dominates 
the pairing at fi/g' 2 S> 1. With screening in matter, we show that in the weak coupling limit the 
Overhauser effect can compete with the BCS effect only for a sufficiently large number of colors. 
Both the BCS and the Overhauser gaps are of order g 4 /fj. in Landau gauge. 
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(N ■ I. INTRODUCTION 

Quantum Chromodynamics (QCD) at high density was first studied in the late seventies [|lj|2j] and has received 
renewed attention, especially in the last three years At large baryon density and small temperature, it is 

I ■ relevant to studies of neutron stars and heavy ion collisions in the baryon rich regime. Theoretical studies of dense 
£SJ ', quark matter have revealed the preferred phase to be a color superconductor with several novel features like color-flavor 
locking, chiral symmetry breaking, color Meissner effect and the formation of a mass gap. 

At asymptotic densities, the Fermi surface for quarks is well defined and the low energy excitations are weakly 
interacting quasiparticles and holes. In QCD, perturbative 1- gluon exchange can provide an arbitrarily weak attractive 
interaction in the color anti-symmetric 3 channel that destabilizes the Fermi surface and causes particle-particle 
pairing. Another mechanism of destabilization is particle-hole pairing (Overhauser effect) studied earlier in a 
variational method and for large N c (|||. More recently, it has been shown that the equations driving the particle- 
hole instability resemble those of the BCS instability, modulo phase-space factors JTof . A novel mechanism for 
destabilizing the Fermi surface through a BCS crystal was also recently suggested in [|11|| and may also compete with 
| the instabilities to be studied below. Its consideration is outside the scope of this work. 

Screening affects the strengths of the BCS and Overhauser pairing. The authors of |l(]] have shown that in the 
regime of weak coupling in 3+ld, the BCS effect is dominant upto a large number of colors. In 1+ld, the drastic 
D . reduction in phase space for BCS pairing with no comparable reduction for Overhauser pairing means that the latter 
' dominates whatever N c . In this work, we wish to establish the case for 2+ld, taking into account the screening in 
, this space-time dimension. 
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This paper is organized as follows. In section II, we derive the quark 1-loop contribution to the gluon propagator 



xi 

in a covariant gauge, at finite temperature and density (T^cC /i for cold and dense matter). Section III motivates 
the symmetry of the particle-particle and particle- hole condensates in 2+1 dimensions, followed by formally differing 
derivations of the gap equations for these two objects. The resulting gap equations are solved in the presence of 
electric screening and Landau damping in dense quark matter. We present our results for the gaps in the Landau 
gauge. The conclusions of our analysis are in section IV. 



II. SCREENING IN 2+1D 

We begin by evaluating the gluon polarization tensor in matter in the imaginary time formalism [ p"2"[ . 

where q , p and p+q denote respectively the 4-momentum of the gluon and the two internal quark lines. At finite tem- 
perature and density, one makes the usual identification ip$ — po = (2n + 1)ittT + /i for fermions and iq^ = q$ = 2niirT 
for bosons (here, n labels the Matsubara frequency). In 2+ld, the Pauli matrices 7^ = 0^ = (173, iS) generate a rep- 
resentation of the Clifford algebra {7^7^} = 2g fiu with = diag(l, —1, —1). In an arbitrary covariant gauge with 
gauge parameter A, the inverse propagator in matter (D -1 )^" and LP" are related as |l3| 
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Evaluating the Matsubara sum using contour integration, omitting the vacuum piece (T and (i - independent), and 
analytically continuing iq^ = uo + ie, we get in the limit T — > (with q =| q |) 
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n L P^, with Pj« = - 214!) and P« = ^ (i,i=l to 3). 



where IT^ = n ± P ± y - 

In deriving Eqs. (^[^), the approximations^ ui, q <C /i have been made since energy transfers in the scattering of qq 
(BCS) and qq (Overhauser) pairs are of the order of the gap, and momentum transfers can be taken to be small for 
forward scattering of quarks. In Euclidean space, the screened gluon propagator in Fcynman gauge reads 



D e ,m(x - y) = 



d 3 q 
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where m E 



m? D refers to the Debye mass and m 2 M to the magnetic scale generated by Landau damping. Their values 
(generalized to Nt flavors) follow from (3) and (4) respectively and are given by m 2 E /g 2 ^L = m 2 D /g 2 [i = Nf/2w and 
'q. Note that the coupling g 2 has mass dimensions. The results are similar to that in 3+ Id for 
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III. CONDENSATES AND GAP EQUATIONS IN 2+1D 
A. Condensates 

Usually, diquarks are taken to condense in the parity even + channel. In 2+ Id, parity is defined by inversion 
about 1 spatial coordinate only, as it is an improper transformation. For the sake of definiteness, we define the parity 
transformation in space as Xi — ► X\, £2 — ► — x-i- The parity operator in spin space is then <j\ which implies that a 
term like qq is not invariant under this transformation. Thus, for the particle-hole pseudoscalar condensate with 2 
flavors, we have the action of the parity operator (P) as 
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dd — > —uu 



dd 
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where c and / denote color and flavor respectively. A scalar condensate may still be written down if we allow parity 

P 

to change isospin as u[d) — > +/ — cr\d(u). Then we have 
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The +/— indicates that there are 4 possible sign combinations we can choose in the combined transformation of the 
u and d quarks. Note that this new definition of parity does not affect the way the pseudoscalar transforms. Now, 
we will see that the relative sign is also determined for us by the particle-particle condensates. 

P 

The charge conjugated spinor transforms under parity as q c — > det(Ap)eri<7 c , where Ap denotes the matrix cor- 
responding to the improper Lorentz transformation P. There is an extra minus sign as compared to the parity 
transformation of the quark spinor. The charge conjugation matrix C is &2 : which is antiymmetric. The attractive 



1 The gaps in 2+ld are power law suppressed, and care is needed in introducing this approximation in our perturbative result. 
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color channel 3 is anti-symmetric and the overall wave function for the BCS pairing must also be anti-symmetric, 
therefore, one is forced to choose a 2 for the flavor matrix. Then, the scalar BCS condensate is 

q c X c 2 a[q = u c d - d c u . (9) 

For this to be even-parity, both u and d should have the same sign in their transformation. On the other hand, a 
pseudoscalar BCS condensate will necessarily have a relative sign between the transformations for u and d. In this 
paper, it is the scalar which is studied. 

These problems arise in odd space-time dimensions due to the lack of a consistent definition for a gamma matrix 
that anti-commutes with all the 7^'s, namely 75. We proceed now to the gap equations for the scalar particle-paticle 
and particle-hole condensates. 



B. BCS Gap equations 

In order to derive the BCS gap equation, we will use the standard Nambu-Gorkov formalism, which introduces 
charge-conjugated spinors to enable the writing of the generating functional of QCD at finite chemical potential in 
a form identical to that at zero chemical potential juj. With charge conjugation defined as ip c = C%p T , the Nambu- 
Gorkov spinor is W = (4>,4> T ) and can be expressed as a A.N c .Nf column vector. Neglecting quark mass effects, the 
form of the gap matrix is 



A#(g) = (X 2 ) ab (a 2 ) l3 C {^-(l + a x q) + A 2 (q)±(l - a X $)} 



(10) 



where A± = ^(1 ± a x q) are the particle and anti-particle projectors respectively in 2+ld. (In the weak coupling 
limit a//i <C 1 with a — g 2 /Ait, the equations of motion allow us to replace a x q by the unit matrix so that only the 
gap Ai remains and the anti-gap A2 drops out. Moreover, the anti-gap is strongly gauge dependent.) The self-energy 
in the Nambu-Gorkov formalism obeys the Schwinger-Dyson (SD) equation 
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Here, S(fc) = S^ 1 (k) — 5 _1 (fc) is the proper self-energy with S(k) as defined in p5[ . T° is the quark-gluon vertex 
which we take to be the leading order result in perturbation theory, 
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We also neglect the wave-function renormalization to leading order in the weak coupling so that the diagonal part of 
the self-energy has only the free quark propagators. Therefore, 

where A = 7oA^7o. Determining 52 1 from Eq. (|T^), inserting the ansatz Eq. ( |l0| ) for the gap, and projecting out the 
positive energy part, we obtain the integral equation 
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\ - \k ■ q)D L (k -q)+ (V(l -k-q)- |K| 2 (1 + k ■ q) j ^ 

where K = q — k,D± = l/(K 2 — m 2 M ),DL = 1/(K 2 — rnjj) and D\ — X/K 4 from the gluon propagator in the 
covariant gauge. We note that only the term for the gap Ai (henceforth referred to as G(q)) is relevant at the Fermi 
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surface and that it's gauge dependence multiplies the factor K° / K 2 which is essentially the amount by which the 
quarks are off mass shell. This is not sufficient to tame the infrared divergence from the 1/K 2 of the gauge-fixing part 
of the propagator in 2+ld, in contrast to 3+ld (see below). An estimate of the relative contributions of the electric, 
magnetic and gauge-fixing terms can now be made. 



1. Gauge fixing contribution {D\) 



It is the momenta on the Fermi surface that make the most contribution to the integral, therefore |q — k| 2 = 
2^i 2 (l — cosO). Setting |qj — A* — Q\\, the integration measure becomes dq 3 [idq\\d9, and the continuation to Euclidean 
space gives 



G(k 3 



Xk I dq 3 ^dq \\G{q 3 )(q 3 - k 3 ) 2 



q 2 + q 2 +G(q 3 ) 2 



d9(l + cos6 
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where n = ((N c + 1) /2N c )(g 2 / (2tt) 3 ) . Performing the angular integration exactly, we find 



G(jfcs) = Aktt J 



dq 3 (j,dq\\G{q 3 ) 
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Integrating over q\ | , the pole from the quasi-particle propagator gives 



G(k 3 ) — Xk^ltt 2 J dq 3 
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\qz-k 3 W{q 3 -k 3 ) 2 + V 
The region of interest is q 3 ,k 3 -C /x, since the quarks are almost on mass shell. In that case, we have 

Aktt 2 f , G(q 3 ) 1 



G(k 3 ) = 



dq 3 



V4 + G(q 3 ) 2 \qs-k 3 \ 
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2. The electric contribution (Dl) 

The angular integration yields 



ktt f dq 3 dq\\G{q 3 ) I yj((q 3 - k 3 ) 2 + m 2 D ) 2 + V((g 3 - k 3 ) 2 + m 2 D ) 



G(k3) 2 M J q 2 + q 2 + G(q 3 ) 2 { ((q 3 k 3 ) 2 + m 2 D ) 1 J ' (W) 

As before, the integration over g|| picks up the pole in the diquark propagator, and the form factor is the physical 
regulator for the energy integral. Note that the integrand shuts off for q 3 >• fi, but this is an unphysical regulation 
since we are far from the Fermi surface physics. For the physical region of q 3 <C ^, we obtain 

G{k 3 )^^ 2 ( - 1 . (20) 

J ^qi + G{q 3 ) 2 ^{q 3 -k 3 ) 2 +m 2 D 

T he value of the Debye mass suggests that it is much larger than the gap, so that we may approximate 
\J (<73 — k 3 ) 2 + m? D = mo, with q 3 and k 3 on mass-shell being of the order of the gap. (We have checked the 
consistency of this approximation with the explicit solution for the gap.) The electric term is expected to be of order 
k/tud- 
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3. The magnetic contribution (D±) 



The magnetic scale is generated by Landau damping through 



» 2 ?0 — kg 

m M = 7 m CT3 7^7 . ( 21 ) 

4 |q-k| 



and the gap equation reads 



r ( k \ = * f dq 3f j,d qil G(q 3 ) r d9(l + cos6>) 

[3) 2j qj+qf + G(q 3 )*J (q 3 -k 3 y + |q-k|2 + \m% ' 1 ' 

As q, k are on the Fermi surface, we may write 

r(k k « /" dq 3 dq n G(q 3 ) [ M rf6>(l + cosfl) 

1 3j 2 ./ + flj + G(, 3 ) 2 7 ((Z3 _ , 3 ) 2 + 2^(1 - cosfl) + >l ^; : -^ Qse) ' 1 j 

The physics of Landau damping guides us in evaluating the angular piece. For forward scattering, 9 is small, in which 
case the angular piece reads as follows 



(q 3 -k 3 ) 2 +v?6 2 +3m 2 D (q 3 -k 3 )/4[i6 ' 

where we have set 1 — cosfl = 2 /2 in the denominator and 1 + cosfl = 2 in the numerator. The denominator contains 3 
pieces, which are, respectively, the gluon energy squared, the gluon momentum squared and the damping term. Since 
the scattering of BCS pairs costs little energy ( typically of the order of the gap ), the first piece may be dropped. 
This is equivalent to saying that the quarks are nearly on mass-shell. The validity of this approximation rests on the 
smallness of the gap which will be checked explicitly. 

The remaining 2 terms give us a bound on the angular region in which magnetic binding of the BCS pairs is 
disturbed by Landau damping. Setting the terms comparable, we find 

_ (fmf,| g 3-M) 1/3 
"min ■ V ^ / 

M 

If 6 < m i n , the Landau damping is active, and is ineffective once 9 > 9 m i n . Energy transfer of the order of the gap 
implies that 9 m i n is small if 

m^ /3 G 1/3 < M • (26) 
This also implies that ignoring (q 3 — k 3 ) 2 in the denominator is valid when 

A<I^ . (27) 

It is magnetic gluon interaction between the quasiparticles that builds up the gap, as in 3+1D |l6| ], so we work in 
the regime 9 > 9 m i n - If we keep only the momentum exchange term in the denominator, we obtain for the angular 
integration 

2 f e — dS_ = 4 

»Je min & 2 6V3 m ^ 3 |g 3 -jfc 3 |i/3 

where fJ-9 ma x — (|q — k|) maa; and we have picked the contribution from the lower limit. The contribution from the 
upper limit is subleading. The dq\\ integration yields 

nth \ 2lTK [ ^Gfaa) 1 (r,^ 

[ 3) ~ 6V3 m V3 J Vqi + G{q3 y \ q3 - fc 3 |i/3 ■ ^) 
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Collecting Eqs. ( |l8| , f20| . p9| ) , the contributions of the electric, magnetic and gauge-fixing terms are respectively 
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The gauge-fixing contribution in covariant gauges is large suggesting that the ladder plus screening expansion of the 
SD equation in aj \x may receive additional contributions due to the the enhanced infrared sensitivity in 2+ Id. A 
reorganization of the expansion maybe needed, with perhaps an emphasis on a physical observable from the outset. 
Below, we will carry estimates of the gap in Landau gauge A = 0. The electric contribution gives an exponentially 
small gap if G < tud- The electric screening implies G ~ yYe~ m - D / K7r , which is exponentially suppressed due to the 
large value of mo/g 2 - (A is the UV cutoff). 

The magnetic term drives the formation of the gap since m 2 ^ 3 G 1 / 3 < mu- The gap equation becomes 



G(fc 3 



dq 3 G(q 3 ) 
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where k* = 2K7r/(6 1 / 3 m 2 ^ 3 ). The integral may be split into two regions, < q 3 < k 3 and k 3 < q 3 < A, where in the 
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Then, Eq. ( |3l| ) may be recast as 
G(q 3 )~ 
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with the scale G(k 3 — Go) = Go- Differentiating Eq. ( p2|) twice with respect to k 3 , we obtain 
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which has the solution 
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Since the Neumann function is complex valued for positive k 3 , C2 — (we are solving for a real gap). Using G(k 3 
Go) = Go and the boundary condition G'(k 3 = Go) = 0, we obtain 
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with Go = (-x) 3 fjjj') ' w here £o = 2.405 is the minimum finite solution to J\{x) = —xJ[(x). Go is therefore the 
maximum BCS gap. Note that the condition Go < mjj is satisfied since 
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The conditions imposed by Eqs. (pq ) and (E7]) are also met because 
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Therefore, the approximations made in evaluating the magnetic and electric pieces are justified. It is easy to check 
that the electric piece is down compared to the magnetic piece by 
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As we work in the Landau gauge, we can repeat the preceding analysis using a simplified propagator 
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and arrive at the same equation for the gap, which in Euclidean space is 
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As expected, the integral over q is dominated by the momenta | q |~ \i and 53 <C u. 



(40) 



C. Overhauser Gap equations 

While it is possible in 3+ld to follow the same formalism for the Overhauser pairing as in the BCS case, the method 
is involved in 2+ld due to the absence of 75, which is crucial in inverting S . Chiral projectors are required to 
pair particles and holes but since chirality cannot be consistently defined in 2+ld, the arguments are more subtle. 
However, there is a straightforward way and that is to use the effective action formalism developed in JToj ] . 

The induced action in Euclidean space is 



d 6 x d 3 y J-( x )D^(x-y)JZ(y)+ / d*x $ d^rf 



(41) 



where <9 M = (<9i, 82,83 + u) and J° is the usual colored current mb^^kf-^. We may now Fierz the J J term into the 
relevant particle-hole channel. The simplified propagator (Eq. (f39|)) is used as we are working in the Landau gauge. 
The Fierzing factors are 



Color: 
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Spin : I 



(42) 



Following p0| , we introduce a hermitian bilocal field S(x, y) to linearize the Fierzed form of the JJ term by using a 
Hubbard-Stratanovich tranformation, as 

exp^a-^) j d 3 x d 3 y Mx)rl>(y)]D(x-v)$(yMx)] 



with 
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The action in the quark fields is now linear and the functional integration may be performed. Before doing so, however, 
we make a simplifying ansatz for the bilocal field 



T,(x, y) — 2 cos 
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a{x — y) = 2 cos 
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d 3 q 
(2tt) 3 



(45) 



where P M = (P^,0) and | Pj? |= 2/i. points in the original direction of one of the quarks. S characterizes a 
standing wave of total momentum 2/i. (The pairing is between a particle and a hole at the opposite edges of the 
Fermi surface.) Introducing fermion fields ^(i-j + l) as independent integration variables as in and performing 
the functional integration over fermions, we obtain 



^(l-J-)^= I* x W°*Z-2? 
8 V N C J V 3 J D(x) J 



d 3 g 
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where Q± — [(i^ + q), (±-j- + (73 — in)]. The determinant is over an (4.N c .Nf) x (4.N c .Nf) matrix. The gap equa- 
tion follows by variation and the result is 
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^ 2 f d 3 q F(q) 



2 J (2ttY { q3 + M_}2 + q 2 + F(g)2 

The measure is d 3 q = dq\\dq±dq3, where q\\ = q P and qj_ = q — <7||P- The form factor F(q) decreases rapidly as we go 
further from the Fermi surface, therefore, dominant contributions to the integral over g|| come from the region where 
q\\ ~ -F(<7||). The gap is approximately constant in the qj_ direction. The q± integration extends to = 
where Fq is the maximum Overhauser gap. Thus, it picks up contributions from a larger region than q|| does. As the 
particles will be placed on mass-shell, q^ and qii are of the same order. The contour integration over q% is performed 

12 ^ O,,, = 9,, f^2 



with the constraint |q±| < 2/ie q = 2/iy^ + F(q\\) 2 . With these observations, the gap equation becomes 

FCfcll^^)^ 3 ^"^ 2 [dq\\-; JM = / A "^(fc ± - g± ,fc||- g ||) . (48) 



v /^+F(Q||) s 



The electric and magnetic terms can be dealt with separately, as in the BCS case. The gap is approximately constant 
in the kj_ direction, so we study the dependence on fc|| alone. 
With electric screening in matter, the gap equation reads 

= ktt [ d % f(g||) = 1 —tan' 1 I = = I , (V.)) 

where k = 3(N C - l)g 2 / '4(2tt) 3 . 

In case of Landau damping of the magnetic mode, the gap equation is 

F(ku) = KTT [ dqn - F{qil) ? ,, C — , (50) 




where c = (61n2 + 2v / 37r)/18. In writing Eq (50), we have expanded in (§™|)|£:|| — <7|| D^/A^ and retained only the 
leading term. This ratio is N c dependent and the expansion turns out to be valid only at sufficiently large N c (see 
Eq (|53|)). Noting that mjj is typically much larger than the gap, the magnetic binding will drive the formation of the 
gap. Proceeding as in the BCS case, we convert the integral equation into a differential one, 

3^F"(fc||) + 4fc||F'(fc||) + K ,^^ = , (51) 
where = (c7TK)/(|rn| ) ) 1 / 3 . Eq. (J5l|) is analogous to Eq (B3J) obtained for the BCS gap. The solution follows 



F 



7/6 ,h(6J±fk- 1/6 ) 



« h) - M»£*s*) C" ' ' (52) 



with F = (^-) 3 (j^j > the maximum Overhauser gap. Note that the condition F < mo is satisfied since 

Fo/m D ~ (a/ixfl 2 . As in the BCS case, Fq/ii ~ (W/i) 2 . These observations justify the neglect of electric screening 
effects. The validity of the approximation in Eq. (p(f) implies that the ratio 



(fmjFo) 1 / 3 _ / 2<kxI \ 1/2 .V 7.24 ^ 1/2 



A± \6c(N c -l)J \{N C -1 



(53) 



should be small. So long as A^ c is sufficiently large, the approximation and the resulting estimate of the Overhauser 
gap can be trusted. If N c does not meet this condition, the binding is disturbed throughout the pairing region by 
Landau damping. With this caveat, we can make a comparison of the maximum BCS and Overhauser gaps (Fq and 
Go). We find 

Jo. ; / v (i-w y f54) 

Go - 2 { iyc (l + l/N c )J ' (54j 
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which shows that in 2+ld, the Overhauser gap is larger than the BCS gap for sufficiently large N c . Finally, an 
assessment of the energy budget for the two competing mechanisms yields 

eov (Aj_Fo)Fo / F \ 2 ( Fq^I" 

(55) 



IV. CONCLUSIONS 



In considering BCS and Overhauser pairing in 2+1 dimensions, we have argued for the form of the condensates 
that we have chosen to study. The gap equations are written down and treated with the effects of electric screening 
and Landau damping in matter. For the case of 2+ld, and to leading order in perturbation theory, we have derived 
modifications to the gluon propagator from screening by quarks, using the imaginary time formalism. The standard 
Nambu-Gorkov formalism is more subtle for the Overhauser case due to the absence of 75 and the effective action 
method is applied instead. A comparison of the BCS and Overhauser gap in Landau gauge, shows that the latter is 
preferred with increasing N c . The energy budgeting shows that BCS pairing is negatively affected by the reduction 
of available phase space for pairing in the lower dimension. The gaps themselves are small, being of order (a 2 ///), 
although the enhanced infrared sensitivity in 2+ld prevents us from a definitive conclusion on their gauge indepen- 
dence. Also, the coupling does not run in 2+ld, so it is not possible to set a perturbative scale, and consequently to 
say when a//x is small. Our perturbative analysis assumes that we can take a/ fi small. It is interesting to note that 
unlike the case of 3+ld or 1+ld, the BCS and Overhauser gaps show power law dependence (g 4 /n) in Landau gauge. 
Therefore, they are more sensitive to higher order corrections from perturbation theory than their corresponding 
forms in 3+ld. A numerical analysis of the gap equations in 2+ld would complement our theoretical approach, and 
would also serve to test the validity of some of the approximations. Finally, it would be interesting to consider the 
effects of the crystalline superconducting phase discussed recently by others P^jnj as a competing instability in this 
space-time dimension. 
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